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ABSTRACT 

We present the Robert C. Byrd Green Bank Telescope discovery of the highly eccentric binary millisecond 
pulsar PSRI1835-3259A in the Fermi Large Area Telescope-detected globular cluster NGC 6652. Timing 
over one orbit yields the pulse period 3.89 ms, orbital period 9.25 d, eccentricity ~ 0.95, and an unusually 
high companion mass of 0.74 M 0 assuming a 1.4 M 0 pulsar. We caution that the lack of data near periastron 
prevents a precise measurement of the eccentricity, and that further timing is necessary to constrain this and 
the other orbital parameters. From tidal considerations, we find that the companion must be a compact object. 

This system likely formed through an exchange encounter in the dense cluster environment. Our initial timing 
results predict the measurements of at least two post-Keplerian parameters with long-term phase-connected 
timing: the rate of periastron advance u> ~ 0.1 0 yr 1 , requiring 1 yr of phase connection; and the Einstein delay 
7gr ~ 10 ms, requiring 2-3 yr of timing. For an orbital inclination i > 50°, a measurement of sin i is also likely. 

PSR J1835-3259A thus provides an opportunity to measure the neutron star mass with high precision; to probe 
the cluster environment; and, depending on the nature of the companion, to investigate the limits of general 
relativity. 

Subject headings: 


1. INTRODUCTION 

Globular clusters (GCs) are efficient producers of low-mass 
X-ray binaries (LMXBs) and their descendant millisecond 
pulsars (MSPs; Papitto et al. 2013, and references therein): 
orders of magnitude more MSPs and LMXBs exist, by mass, 
in clusters than in the Galactic field (Camilo & Rasio 2005). 
The dense GC environment heightens the probability of stel¬ 
lar interactions (parameterized by encounter rate y; Verbunt 
& Freire 2014), increasing the likelihood of forming new bi¬ 
naries and of existing binaries gaining new companions. Sys¬ 
tems that rarely (if ever) form through known binary evolu¬ 
tionary channels in the field can in principle form through 
such stellar interactions in GCs, for example: sub-ms pul¬ 
sars; highly eccentric binaries; or unusual binary systems 
like MSP-main sequence (Pallanca et al. 2010, and references 
therein), MSP-MSP, or MSP-black hole (MSP-BH) binaries 
(Ransom 2008). Such systems would allow astrophysical 
studies that may not otherwise be possible, e.g., strong-field 
tests of gravity with MSP-MSP or MSP-BH binaries. 

The Fermi Large Area Telescope (LAT) has found MSPs 
to be nearly ubiquitous y-ray emitters; GeV emission from 
GCs (Abdo et al. 2010; Tam et al. 2011) may originate from 


the clusters’ MSP populations (e.g.. Venter et al. 2009). The 
LAT-detected GCs included NGC 6388 and NGC 6652 (Abdo 
et al. 2010), neither of which contained any known MSPs, 
but whose y-ray luminosities implied large MSP populations. 
NGC 6388 is particularly interesting due to its high encounter 
rate (e.g., Freire et al. 2008; Maxwell et al. 2012); NGC 6652 
may also have a higher encounter rate than previously thought 
(Noyola & Gebhardt 2006). The presence of a MSP popula¬ 
tion is supported by the number of X-ray sources, including 
LMXBs, in both clusters (at least two in NGC 6652 and five 
in NGC 6388: Maxwell et al. 2012; Stacey et al. 2012). We 
searched these clusters for radio pulsars; here we report on 
the discovery and timing of a highly eccentric binary MSP in 
NGC 6652. 

2. OBSERVATIONS AND PULSAR SEARCH ANALYSIS 

We observed NGC 6388 and NGC 6652 (Table 1) with 
the National Radio Astronomy Observatory 1 Robert C. Byrd 
Green Bank Telescope (GBT) and the Green Bank Ultimate 

1 The National Radio Astronomy Observatory is a facility of the National 
Science Foundation operated under cooperative agreement by Associated 
Universities, Inc. 
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Pulsar Processing Instrument (GUPPI 2 ) backend (DuPlain 
et al. 2008), at S band (2 GHz) in coherent search mode. 
The data were taken with 2048 spectroscopic channels and 
an effective bandwidth of 700 MHz (accounting for excised 
radio frequency interference), with coherent dedispersion at 
the clusters’ predicted dispersion measure (DM) values (Ta¬ 
ble 1). The observing log and minimum detectable flux densi¬ 
ties S m in for an assumed 10% duty cycle are given in Table 2. 

The data were analyzed using presto 3 (Ransom 2001). 
Time series were dedispersed at 1900 DMs from 0- 
691 cm 3 pc for NGC6652 and 5456 DMs between 0- 
800 cm 3 pc for NGC6388, and were searched for periodic¬ 
ities. We searched for accelerated signals over z = ±200 
Fourier bins (cf. Ransom et al. 2001), corresponding to max¬ 
imum line-of-sight accelerations between + 300-3000 cm s 2 
for a 5 ms pulsar. 

3. DISCOVERY AND INITIAL TIMING ANALYSIS OL PSR 
J1835-3259A 

We discovered PSR J1835-3259A (hereafter NGC6652A) 
in the direction of NGC6652 (DeCesar et al. 2011), with 
the fundamental frequency at an acceleration of 11.1 cm s 2 
(z = 9). Figure 1 shows the presto discovery plot, and Ta¬ 
ble 2 contains estimates of the 2 GHz flux density S 2 . We 
discuss the unexpectedly low DM value of 63.35 cm~ 3 pc be¬ 
low (Section 4.1). 

We fit the Doppler-shifted P and P (Table 2) with a phase- 
incoherent orbital model Freire et al. (2001), using a routine 
by R. Lynch (private communication) employing mpfit 4 , and 
found a very eccentric orbit (e > 0.7). Starting with this or¬ 
bital model, we ran tempo 5 iteratively on the pulse times of 
arrival (TOAs; Table 2) to converge on a family of timing so¬ 
lutions. We phase-connected the first five observations; we 
did not observe the pulsar at periastron, between observations 
5 and 6, so we allowed the phase between these observations 
to remain arbitrary (i.e., we kept a JUMP between these ob¬ 
servations’ sets of TOAs). Using the DD model (Damour & 
Deruelle 1985, 1986), we find e = 0.968. “Faking” phase 
connection by removing the JUMP yields e = 0.950; alterna¬ 
tively, including arbitrary phase JUMPs between all TOA sets 
yields e ~ 0.8, which we take to be the lowest possible e of 
this system. 

The best-fit DD timing model parameters are in Table 3, with 
fit residuals in Figure 2. The systematics in the residuals are 
present in all our fits, including those with JUMPs between all 
observations; we attribute them to parameter, and therefore 
phase, uncertainties resulting from the lack of TOAs through 
periastron. We stress that the initial timing parameters in Ta¬ 
ble 3 belong to a family of solutions—a unique determination 
of the MSP’s timing solution requires further observations, 
especially through periastron passage. If further timing con¬ 
firms the parameters, then NGC6652A will be the most ec¬ 
centric binary MSP known to date. 

4. DISCUSSION 

We adopt timing parameters from the e = 0.950 model, 
pulsar mass m p = 1.4 M 0 , and cluster parameters from Har¬ 
ris (1996, 2010 edition) for all calculations, unless otherwise 
stated. 

2 http://dx.doi.org/10.1117/12.790003 

3 http://www.cv.nrao.edu/ sransom/presto/ 

4 http://cars9.uchicago.edu/software/python/mpfit.html 

5 http://tempo.sourceforge.net/ 


4.1. Cluster Membership 

The discrepancy between the discovery and predicted DMs 
(63.35 and 190cm 3 pc, respectively) initially led us to ques¬ 
tion the MSP’s cluster association (DeCesar et al. 2011). 
However, the Cordes & Lazio model commonly has uncer¬ 
tainties of a factor 0.5-2, and sometimes larger. The mea¬ 
sured DM is consistent with the low optical reddening Lg-v = 
0.10 + 0.02 (Ortolani et al. 1994) and estimated X-ray absorp¬ 
tion column Nu ~ 5.5 X 10 20 cm 2 (Predehl & Schmitt 1995) 
toward NGC 6652. 

The high e of NGC6652A is much more probable in a GC 
than in the field due to the high probability of stellar encoun¬ 
ters (Camilo & Rasio 2005), discussed further below. Addi¬ 
tionally, given the beamwidth of the GBT at S-band (6.3'), we 
estimate a ~ 0.2% chance of finding an unassociated MSP co¬ 
incident with NGC 6652 (assuming an isotropic distribution 
of known galactic MSPs). We conclude that the MSP is al¬ 
most certainly a cluster member. 

4.2. Nature of the Companion 

The minimum companion mass (orbital inclination i - 90°) 
is m Ci min — 0.74 M 0 (Table 3). Comparing with the Australia 
Telescope National Facility (ATNF) Pulsar Catalog 6 shows 
that the companion is unusually massive; it may be a main 
sequence (MS) or evolved star, or a compact object. Based on 
the cluster’s age (11.7 ± 1.6 Gyr; Chaboyer et al. 2000), the 
main-sequence turn-off mass is « 0.8 M 0 (Stacey et al. 2012). 
For i < 70°, m c > 0.8 M 0 , limiting the range of inclinations 
for which an unevolved MS companion is possible (c.f. Freire 
et al. 2007). 

For a non-compact companion, significant tides at perias¬ 
tron will circularize the orbit. The circularization, or dissipa¬ 
tion, timescale to for an eccentric binary system can be esti¬ 
mated as (Socrates et al. 2012, and references therein) 

m c a p 

= 3k L rGw 2 R 5 

where ap = a sin i (1 - e 2 ), r is the constant tidal lag time of 
the companion, k[ is the Love number, and R is the compan¬ 
ion’s radius. For high e, the tidal quality factor Q is related 
to t by Equation 23 of Socrates et al. (2012). We estimate k\ 
to be between 0.05-0.15 for both MS and WD companions, 
based on calculations with Modules for Experiments in Stel¬ 
lar Astrophysics (mesa; Paxton et al. 2011; Brooker & Olle 
1955). For Q = 10 6 , the circularization timescales are ~ Myr 
for a MS companion and ~ 10 4 Gyr for a WD companion. We 
conclude that the companion is a compact object, whose mass 
and nature will be constrained through further timing. 

The merger timescale from gravitational wave-driven inspi¬ 
ral depends on e as /merge (1 - e 2 ) 1 ^ 2 (Peters 1964). For 
i = 90 (iHc.min = 0.74 M 0 ) and e = 0.95, / mer ge ~ 12 Gyr; vary¬ 
ing e yields a range of / mer ge ~ 1 Gyr (e - 0.975) to > 100 Gyr 
(e < 0.9). The system may therefore be disrupted (Sec¬ 
tion 4.4.2) before it has time to merge. In the event of a 
merger, when the system comes into contact, the outcome will 
depend on the exact nature of the binary. Stable mass trans¬ 
fer will be possible for q = m c /m p < 2/3 (i > 50°), form¬ 
ing an ultra-compact X-ray binary and possibly an isolated 
MSP. For larger inclinations, the mass transfer will be unsta¬ 
ble; while a black hole would form from accretion-induced 
collapse (AIC; Giacomazzo & Perna 2012, and references 

6 http://www.atnf.csiro.au/people/pulsar/psrcat/ 
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therein) if the system mass exceeds the maximum NS mass, 
it is unclear whether substantial mass would be ejected from 
the system during unstable mass transfer, preventing AIC (L. 
Bildsten, private comm.). An eventual merger of this system 
may result in a long-GRB-like, calcium-rich transient (e.g.. 
King et al. 2007) if the companion is a massive WD, or a 
short GRB for a NS companion (e.g., Grindlay et al. 2006). 

4.3. Post-Keplerian Parameters and Mass Constraints 

Finding m c and m v requires measurements of at least two 
post-Keplerian (PK) parameters (we employ the general rel¬ 
ativistic formalism of Damour & Taylor 1992). Our prelim¬ 
inary timing solution predicts that the rate of change of the 
longitude of periastron passage to (i.e., the orbital precession 
rate) is to > 0.08° yr ', the Einstein delay yc,R > 10 ms, and 
Pi, > 6x 10~ 12 s s -1 . Because of the high e , we will measure to 
with high precision: from simulations assuming the e = 0.95 
orbital model parameters, we find that we will measure to with 
> lOOcr significance after one year of timing, yielding the to¬ 
tal system mass ( to oc M^ 3 , where M tot = m p + m c ) and con¬ 
straints on OT p and m c . Knowledge of the pulsar position (re¬ 
quiring 1 yr of timing or an interferometric detection) would 
yield a measurement of to with one month of phase-connected 
timing. Our simulations also show that ycR will be measured 
with 10% uncertainty with 2.5 yr of phase connection. 

We may also measure one Shapiro delay parameter, s = 
sin;. For i > 50°, the timing residuals from s are signifi¬ 
cantly larger than the ~ 20 ps uncertainties in the pulse TOAs 
we used to build the timing model. Statistically, it is most 
likely that the MSP companion is a WD, requiring i > 40° 
for ;h c < 1.4M 0 . Even a marginal detection of Shapiro delay 
will yield a precise s because to and .v are nearly orthogonal 
in the mass-mass diagram (Lynch et al. 2012). With these 
two PK measurements, we would precisely measure m p and 
m c . We note that the very precise mass of PSR J1807-2500B 
(NGC6544B) was measured in this way (Lynch et al. 2012). 

4.4. System Origin 

The vast majority of field MSP binaries have circular or¬ 
bits (cf. Champion et al. 2008) from dissipation during the 
mass-transfer phase (Phinney 1992); known eccentric systems 
in the field are either double NSs (with eccentricity coming 
from a second supernova (SN) kick; e.g., Brandt & Podsi- 
adlowski 1995), disrupted triples (Champion et al. 2008), or 
possibly NS-He WD binaries with circumbinary disks (Anto- 
niadis 2014, and references therein). In contrast, a number 
of the MSP binaries in GCs are substantially eccentric (e.g., 
Freire et al. 2008; Lynch et al. 2012) 7 , with a likely origin 
in dynamical encounters (e.g., Verbunt & Freire 2014). The 
highest-e binary MSP currently known, PSR J0514-4002A 
(NGC1851A; Freire et al. 2007), has e = 0.888 and an 
unusually massive (m c > 0.96 M 0 ) companion, similar to 
NGC6652A. Here we consider the plausibility of several 
mechanisms through which NGC6652A could have gained its 
high e. 

4.4.1. Possible Formation Mechanisms 

There are several ways to form a high-eccentricity system 
like NGC6652A. An initially circular orbit may gain eccen¬ 
tricity from 3-body encounters with other stars in the GC (Ra- 
sio & Heggie 1995). For a double neutron star (DNS), the 

7 Also see http://www.naic.edu/-pfreire/GCpsr.html. 


eccentricity could have been imparted on the system by the 
SN kick of a massive companion (Brandt & Podsiadlowski 
1995, and references therein). The system could also have 
formed through an exchange encounter, in which the original 
companion was ejected from the system and the third body 
became the new companion (e.g., Verbunt & Freire 2014). In 
this case, the new companion can be any type of compact ob¬ 
ject. 

The SN kick is ruled out by observational evidence that all 
known radio pulsars with NS companions have spin periods 
of 20-100 ms (Tauris 2011), suggesting that MSPs cannot be 
fully recycled by short-lived, massive companions. The first 
scenario is plausible, as using Equation 5 of Rasio & Heggie 
(1995), we find that » 11.8 Gyr (comparable to the GC age) of 
non-exchange three-body interactions would be needed for a 
binary in an initially circular orbit to gain e = 0.95. However, 
the exchange encounter scenario seems most natural, and we 
discuss this mechanism in more detail below. 

Other scenarios for the origin of the binary’s eccentricity 
include a physical collision between a MSP and a giant star 
(e.g., Freire et al. 2007), or a triple system in which the out¬ 
ermost companion is pumping the eccentricity of the inner 
binary (e.g., B1620-26; Thorsett et al. 1999). These mecha¬ 
nisms cannot be excluded a priori , but are outside the scope 
of this letter. 

4.4.2. Dynamical Formation Through an Exchange Encounter 

We consider a dynamical encounter resulting in a compan¬ 
ion exchange, using fewbody (Fregeau et al. 2004) to simu¬ 
late a particular scenario. As a progenitor system, we take 
the current most common MSP binary in GCs: a MSP in 
a circular 2d orbit with a low-mass companion. We chose 
m c — 0.2 Mq, which follows from a binary period of 2 d using 
the period-core mass relation from Tauris & Savonije (1999) 
for Pop II stars. We simulated 5000 encounters between this 
binary and a third body, drawing the incoming velocities from 
a Maxwellian distribution 8 using cr = lOkms -1 distributed 
between 0 and 30a (where a is the binary’s semi-major axis). 
For the third body we assume a WD with 0.7 M 0 . 

Approximately 70% of the encounters result in an ex¬ 
change, with the low-mass companion ejected and an eccen¬ 
tric binary remaining. The new binary has a range of eccen- 
tricies strongly biased toward high values, with 66% of the 
new systems having e > 0.8, but energies comparable to that 
of the progenitor. The orbit has expanded due to the factor of 
3.1 increase in m c , leading to a factor of 3.1 7 2 increase in 1\. 
Therefore, systems with e w 1 and P\, ~ 10 d are naturally 
formed through this mechanism. If no exchange happened, 
then the binary remains close to the 2 d initial period, albeit 
with enhanced eccentricity. 

We estimate the frequency of encounters between a particu¬ 
lar NGC6652A-like binary and a single star in NGC 6652 us¬ 
ing the single-binary encounter rate y from Verbunt & Freire 
(2014), normalized to M4, and find y 6652 ~ 6.7yM4- A 
NGC6652A-like binary in M4 would encounter single stars 
at a rate £ 1+2 ~ 0° C ,M4/To)c r i+2VM4, where p c is the GC core 
density and oq +2 is the gravitationally-focused single-binary 
cross-section (Equation A2, Leigh & Sills 2011). In M4, 
this encounter rate is £i+ 2 ,m 4 ~ 0.17 Gyr 1 ; in NGC 6652, 

8 While the velocity profile of NGC 6652 has not been measured directly, 
we estimate a velocity dispersion cr of about 10-15 km s _1 , scaling from glob¬ 
ular clusters with similar physical core radii (NGC 6388, 6093, and 6441); 
comparable values are obtained by McLaughlin & van der Marel (2005). 
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i'] + 2,6652 ~ 1 Gyi' 1 . We note that the core radius of 1".15 
measured by Noyola & Gebhardt (2006) is much smaller than 
that from Harris (1996, 2010 edition), yielding 76652 ~ 38yM4 
and f 1 + 2,6652 ~ 6 Gyr 1 . The companion exchange scenario is 
therefore quite plausible. The position of NGC6652A in the 
cluster may give additional clues to its formation (cf. Phinney 
& Sigurdsson 1991). 

5. CONCLUSIONS 

We discovered one new MSP, NGC6652A. Although 
NGC6388 and NGC6652 are expected to host substan¬ 
tial MSP populations, cluster MSPs are extremely faint— 
detecting them requires long integration times and the largest 
telescopes in the world. We did not find more MSPs in these 
GCs simply because we are sensitivity limited. NGC6652A 
is an intriguing source for southern-hemisphere Square Kilo¬ 
meter Array (SKA) precursors and eventually the SKA Mid- 
Frequency Aperture Array 9 . 

Our timing analysis over 1.2 orbit of NGC6652A shows 
that the MSP is in an extremely eccentric binary system with 
an unusually massive compact companion. The system quite 
plausibly formed through an exchange encounter in the dense 
cluster environment. We cannot exclude all other formation 
mechanisms a priori ; determining the nature of the compan- 

9 https://www.skatelescope.org/mfaa 


ion will help discriminate between scenarios. Similarly, a pre¬ 
cise position will help determine whether the binary is dynam¬ 
ically relaxed (and hence close to the core as expected from 
mass segregation) or has been kicked out of the core by a re¬ 
cent encounter. 

With a phase-connected timing solution over > 1 yr, we will 
uniquely determine the MSP’s timing parameters and measure 
its position and co to high precision. After ~ 2.5 yr of timing, 
we expect to measure ycR, allowing measurements of m p and 
m c and clarifying the companion’s nature. If i > 50°, we may 
also measure sink New timing observations are underway, 
and will be reported upon in a subsequent paper. 
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TABLE 1 

Targeted Globular Clusters" 


Cluster Name 

a 

5 

1 

b 

Distance 

Predicted DM 




(degrees) 

(degrees) 

(pc) 

(cm -3 pc) 

NGC6388 

17 h 36 m 17“88 

-44° 44'O'.'24 

345.56 

- 6.74 

11.6 ± 2.0 fc 

340 r 

NGC6652 

18 h 35 m 44)86 

-32°59'25'.'10 

1.53 

- 11.37 

9± \ d 

190 r 


" The cluster positions were set to the optically determined positions of the cluster centers. 
b Moretti et al. (2009) 
c Cordes & Lazio (2002) 
d Chaboyer et al. (2000) 


TABLE 2 
Observation Log" 


Date 

tint 

(s) 

s m J 

(/'Jy) 

S/N 

S 2 C 

Wy) 

n toa 

f TOA 

(s) 

(°toa) 

(its) 

B ary centric P 
(ms) 

B ary centric P 
(10 -12 ss -1 ) 

z 

(Fourier bins) 

NGC6652 

2010 Oct 19 

9470 

5.8 

25.2 

22.5 

8 

1200 

20 

3.88937447(3) 

-1.223 ±0.027 

8.50 

2010 Oct 21 

10062 

5.6 

14.1 

12.2 

8 

1250 

22 

3.88915225(3) 

-0.888 ± 0.023 

7.00 

2010 Oct 22 

8878 

6.0 

5.8 

5.4 

7 

1100 

17 

3.88904917(3) 

-0.916 ±0.028 

6.00 

2010 Oct 21 d 

7504 

6.5 

17.4 

17 

4 

940 

23 

3.888933(1) 

-1.33 ±2.09 

. . .(I 

2010 Oct 24 

6701 

6.8 

14.2 

15.1 

8 

840 

21 

3.88878963(5) 

-1.883 ±0.060 

6.25 

2010 Oct 29 

3086 

10.1 

12.6 

19.7 

7 

380 

32 

3.8891541(1) 

-0.934 ± 0.286 

0.50 

NGC6388 

2010 Oct 21 

6278 

7.3 









2010 Oct 24 

3403 

9.9 









2010 Oct 29 

4080 

9.0 









2011 Feb 05 

5412 

7.8 









2011 Apr 08 

6130 

7.4 









2011 May 06 

4338 

8.7 










" All observations were taken at 2 GHz with ~ 700 MHz effective bandwidth and 40.96 /is time resolution. The beamsize was ~ 6'. 
b ‘S’min was calculated with the radiometer equation for pulsed signals (Appendix A 1.4 of Lorimer & Kramer 2005), using S/N m j n = 
5 and pulse width W = 0.1 P. For observations at S band, the correction factor = 1.05, gain G = 1.9KJ -1 , and T rec ^ 22 K; for 
these observations, the number of polarizations n p = 2, and effective bandwidth A/ = 700 MHz. The sky temperature T s k y = 0.8 K 
for NGC6652 and 1.4 K for NGC6388 (de Oliveira-Costa et al. 2008, 2010), giving T sys = T rec + T, ky = 22.8 K and 23.4 K, 
respectively. 

c The 2 GHz flux density S 2 was calculated with the same radiometer equation parameters as for S m j n , but using the measured S/N 
rather than S/N m j n = 5 and the measured W = 0.0625P rather than W = 0.1P. The smaller measured W explains why we find 
S 2 < S m in on 2010 October 22. 

d The DM used for coherent dedispersion of individual channels was accidentally set to 9.0 cm -3 pc for this observation. As a 
result, the MSP was not found in a straightforward acceleration search. 
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TABLE 3 

NGC6652A Timing Solution" 


Timing Parameter 

With JUMP 

Without JUMP 

Right Ascension* (J2000.0) 

18 h 35 m 44 s .856 

18 h 35 m 44 s .856 

Declination* (J2000.0) 

-32° 59' 25".08 

-32° 59' 25" .08 

Dispersion Measure 0 (cm -3 pc) 

63.35 

63.35 

Spin period, P (ms) 

3.888824(1) 

3.8888289774(4) 

Spin period epoch (MJD) 

55488.931354 

55488.931354 

Spindown rate rf , P (s s -1 ) 

0 

0 

Orbital period, Pb (days) 

9.2460(5) 

9.2459(5) 

Projected semimajor axis, x (s) 

19.6(3) 

19.09(5) 

Eccentricity, e 

0.968(5) 

0.950(1) 

Epoch of periastron passage. To (MJD) 

55477.061(5) 

55477.0401(6) 

Longitude of periastron, to (degrees) 

291(1) 

289.2(2) 

Minimum companion mass, m cxnm (M 0 ) 

0.765(14) 

0.736(3) 

Fit x 2 per degrees of freedom 

347.9/35 

369.37/36 


" The solution uses the DD model (Damour & Deruelle 1985, 1986) and the TDB 
time system. The error on the last digit(s) of each parameter value is denoted 
in parentheses. The true solution is one in a family of solutions represented 
by the parameters listed here. The middle column gives the timing parameters 
obtained when allowing an arbitrary number of pulsar rotations between the 
fifth and sixth observations. The right column gives the parameters obtained 
with “forced” phase connection. The unique solution will be determined with 
further timing observations that include a periastron passage. 

* The position was fixed at the cluster’s center. 

‘ The DM was fixed to this best value from the discovery observation (2010 
October 19). 

d The spindown rate was fixed at zero; a phase-connected timing solution span¬ 
ning ~ one year will measure this parameter. 


2 Pulses of Best Profile 


Candidate: 
Telescope: 
Epoch to 
Epoch b 


ACCEL_Cand_1 

GBT 

, = 55488.93135416666 
, = 55488.93029454615 


Search Information 

ra J2000 = 18:35:44.8560 DEC J2000 

Best Fit Parameters 

Reduced X* = 3.483 P(Noise) < 3.59e—1 
Dispersion Meosure (DM) = 63.350 

P, 


= -32:59:25.0800 


topo ' 





(ms) = 3.889374474(33) P„ (ms) - 
P 'VoV= (‘/’l- -1 -223(27)*IT' 2 P’*, (s/s) - 
P"«p 0 6/s 2 > - 0.0(1,8)«10 (s/s 2 ) ■ 

Binary Parameters 

P orb ( S ) = N / A 


(*8.9a) 

3.889014266(33) 
-1.499(27)x10 - 2 
« 0.8(1.8)x10 17 



3 2 10 
Reduced g 2 


5x10"® 0 -5x10“® 

Period - 3.88937457 (ms) 


Fig. 1. — presto discovery plot. Top left: Two cycles of the summed pulse profile. Left: Two pulse cycles, with the signal in 64 sub-integrations split evenly 
over the full integration time. The x 2 plot on the right shows the steadily increasing S/N over time. Middle: Two pulse cycles with the signal in 32 frequency 
sub-bands split evenly over the full bandwidth. Bottom middle: The x~ increases dramatically near the pulsar DM. Right: The maximum^ 2 determines the best 
P (top) and P (middle). Bottom right: Covariance between P and P. 
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Fig. 2. — Top: Predicted spin periods for a circular orbit (dashed line) and the DD orbital model with a JUMP between observations 5 and 6 (e = 0.967; solid 
line). The incoherently measured P and P, including uncertainties, are overlaid in black. Middle: Residuals from subtracting P predicted by the e = 0.967 model 
from the measured P. Bottom: Timing residuals from TOAs fit with the e = 0.967 model. We attribute the systematics in the TOAs to phase uncertainties which 
will be resolved with a fully phase-connected timing solution. 
















